Multi-tiered quantization of channel state information in multiple antenna systems

ABSTRACT

A multi-tiered CSI vector quantizer (VQ) is provided for time-correlated channels. The VQ operates by quantizing channel state information by reference to both the current channel state information and a prior channel state quantization. A system is also provided that uses multi-tiered CSI quantizers. Enhanced signaling between the transmitter and receivers is provided in order to facilitate the use of multi-tiered CSI quantizers.

BACKGROUND

One of the most promising solutions for increased spectral efficiency in high capacity wireless systems is the use of multiple antennas on fading channels. The fundamental issue in such systems is the availability of the channel state information (CSI) at transmitters and receivers. In general, if the receivers and transmitter have an access to CSI, the system throughput can be significantly increased. While it is usually assumed that perfect CSI is available at the receivers, the transmitter may only have partial CSI available due to the feedback delay and noise, channel estimation errors and limited feedback bandwidth, which forces CSI to be quantized at the receiver to minimize feedback rate. There is described here an improvement in the quantization of channel state information in a multiple antenna system.

SUMMARY

A multi-tiered CSI vector quantizer (VQ) is provided for time-correlated channels. The VQ operates for example by quantizing channel state information by reference to both current channel state information and a prior channel state quantization. A system is also provided that uses multi-tiered CSI quantizers. Enhanced signaling between the transmitter and receivers is provided in order to facilitate the use of multi-tiered CSI quantizers. These and other aspects of the device and method are set out in the claims, which are incorporated here by reference.

BRIEF DESCRIPTION OF THE FIGURES

Embodiments will now be described with reference to the figures, in which like reference characters denote like elements, by way of example, and in which:

FIG. 1 is an illustration of a channel vector space according to known principles in the art;

FIG. 2 is an illustration of the channel vector space of FIG. 1 with a more fine-grained quantization;

FIG. 3 is an illustration of the channel vector space of FIG. 2 with a multi-tier quantization;

FIG. 4 shows the structure of a quantization system;

FIG. 5 shows the operation of an algorithm for designing a multi-tiered quantizer;

FIG. 6 shows a region of one codebook and the corresponding regions of the corresponding next higher tier codebook for use in a multi-tiered quantizer;

FIG. 7 shows an example of the operation of a multi-tiered quantizer;

FIG. 8 is a flow diagram showing the operation of a multi-tiered quantizer;

FIG. 9 shows the operation of the quantizer;

FIG. 10 shows typical eigenmode and singular value coherence times; and

FIG. 11 shows the results of a simulation comparing a multi-tiered quantizer with non-tiered quantizers.

DETAILED DESCRIPTION

In a multiple antenna system as for example shown in FIG. 4, information is transmitted over multiple channels 34 corresponding to multiple antennas 36. Each channel 34 has a state that affects the propagation of information over the channel. The state of multiple channels 34 between a transmitter 32 and one or more receivers 30 in a multiple antenna system can be expressed as a vector. As the channel state changes, this vector moves through the channel vector space. The channel vector space may be separated into regions (see FIG. 1). Each region may be represented by an index. If each region corresponds to the part of the space closest, by some metric, to a particular member of a set of points in the space, then the regions are known as Voronoi regions 20 and the points are known as centroids 22. In order to maximize throughput, it is preferred to associate each index to a centroid 22, which represents the Voronoi region 20 which is the part of space closer to that centroid 22 than any other.

A vector quantizer (VQ) with multi-tiered quantization is aimed at transmission channels with memory, in which there is a need to reduce the feedback bandwidth and allow the system to automatically adjust the quantizer resolution to the rate of channel changes. In an exemplary design of a multi-tiered CSI vector quantizer (i.e., the description of centroids and Voronoi regions) for multiple-input, multiple-output (MIMO) channels with memory, the VQ uses multiple optimization steps.

In the typical CSI VQ, the quantization of the channel vector space can be illustrated as in FIG. 1: the CSI space is tessellated by Voronoi regions 20 with corresponding centroids 22 that represent all vector realizations within each Voronoi region (for each centroid, the corresponding Voronoi region is the set of points closer to that centroid than any other, according to some metric). The number of such regions (centroids) is defined by the number of available bits, which also influence the quantization error of the VQ. FIG. 1 shows the situation when the channel CSI is correlated in time and follows some trajectory 28 in time.

The quantization error can be decreased if the CSI VQ resolution is increased using more bits in feedback link. FIG. 2 shows the identical trajectory 28 of the channel vector realization as in FIG. 1 with increased number of centroids 22 and Voronoi regions 20A. In FIG. 1 only two different centroid indices would be used to represent the channel trajectory, whereas in FIG. 2, four such indices would be used, which would result in more precise representation of the actual channel changes. The price for such improvement is a larger number of bits needed to characterize the quantized CSI indices that must be fed back to the transmitter.

The typical trajectory 28 of CSI vectors is partially predictable in a sense that the channel realizations between consecutive transmission epochs within the same frequency band are correlated. The correlation increases with decreasing relative speeds of receiver-transmitter pairs with the net effect of trajectories being statistically contained within a given Voronoi region for a predictable amount of time. The time metrics may be quantitative described in various ways such as using time metrics called eigenmode coherence time and singular value coherence time. In CSI VQ context, the longer the coherence time, the less frequent the changes in VQ indices that need to be reported back to the transmitter.

A multi-tiered VQ allows for a significant reduction of the feedback rate for systems in which channel coherence times are fairly long. During the design of the quantizer, the Voronoi regions are optimized according to any chosen criterion in 2, 3, 4 and more tiers, in which consecutive Voronoi regions are embedded in the previous ones as shown for example in FIG. 3 for a 2-tiered design.

In the example of FIG. 3, a 2-tiered CSI VQ is divided into primary 23 and secondary 24 Voronoi regions and corresponding centroids 25, 26. In the first phase of the VQ operation, only the primary regions are used to assign the primary centroid indices to the channel vectors. In the second phase, only the secondary centroids and the primary centroid within the first identified primary region are reported to the base station until the channel vector realization leaves the primary region. In this way, as long as the channel vector does not change very rapidly, high quantization resolution can be obtained at much lower feedback rate as the quantization points are concentrated within a space of single primary Voronoi region. Moreover, this mechanism allows the receiver to automatically adjust the vector resolution to the rate of channel changes. The transmitter and receiver must have a way of identifying for which Voronoi regions (primary, secondary etc.) the VQ indices are reported.

The following notation is used in describing an exemplary multi-tiered VQ:

-   -   M—the number of tiers in the CSI vector quantizer design     -   m_(k)—the current tier index at receiver k     -   m^(k)—the current base station tier index of receiver k     -   N_(m)—the number of bits for CSI representation at each tier of         the CSI MIMO VQ.

A system using a dual VQ codebook design for quantization of channel state information in a multiple antenna system is shown in FIG. 4. This example shows a system according to the inventors' U.S. patent application Ser. No. 11/754,965 filed May 29, 2007. A multi-tiered VQ may be used for eigenmode and singular value codebooks in systems ranging from only one active receiver at a time to systems with multiple receivers being active simultaneously (where we define being active as receiving transmissions). The design of the multi-tiered codebooks can be applied to matrices of orthogonal eigenmodes, subsets of eigenmodes and scalar singular values as necessary. The following descriptions may be applied to any type of CSI quantizing solution.

In FIG. 4, a transmitter 32 communicates with a receiver 30 over a feedforward channel 66 and a feedback channel 38 using antennas 36. The receiver 30 includes a channel estimator 40, linear processor 68 for decoding a transmission, a singular value processing unit 42, a power allocation and eigenmode selector 44, and codebooks 48 and 46. The receiver 30 may use various known electronic processors for its parts, and in one embodiment may use a monolithic application specific chip. The functions of the receiver 30 may be provided partly or entirely by hardware, firmware and/or software. The transmitter 32 includes an indexer and optimizer 54 and stored modulation and power allocation matrices 56 and 58 respectively. An input data stream 60 is fed to a modulator 62 that applies a linear modulation matrix selected from the stored modulation matrices 56. The modulated data stream is fed to a power allocator 64, which applies a power allocation matrix selected from the stored power allocation matrices 58. The system of FIG. 4 works as follows:

-   -   1. Before the transmission epoch, each receiver 30 estimates 40         its channel matrix H 34 for the feedforward channel 66 and uses         this information to perform 42 the singular value decomposition         (SVD) of the matrix.     -   2. The eigenmode and singular value (power allocation)         components are separately quantized 44 using two codebooks V 46         and D 48, respectively.     -   3. The indices 50 of the selected codewords are fed back to the         transmitter using a feedback channel 38.     -   4. The transmitter uses all the indices from all receivers 50,         52 in the system to choose 56, 58 the pre-computed linear         modulation and power allocation matrices B 62, 64 and S,         respectively. The choice is based on a predefined set of rules         (maximum throughput, fairness, etc.).     -   5. The signal (x₁ to x_(N) _(T) ) 60 is modulated using the         selected linear modulation and power allocation matrices B 62         and S 64 and transmitted via the feedforward channel 66.     -   6. The transmitted modulated signal is processed by the receiver         68.         The transmitter 32 thus has a processor configured to carry out         the above steps 4-5 and each receiver has one or more antennas         36 and a processor configured to carry out steps 1-3 and 6.

Referring to FIG. 5, design of the multi-tiered codebooks (D or V) is performed as follows:

-   -   1. Based on the desired system parameters (types of channels,         required feedback rate required performance etc.) set parameters         M and N_(m), for each value of m=1, 2, . . . , M.     -   2. Set m=1 (step 70).     -   3. Any of various vector quantizers may be used to design a         receiver VQ using N₁ resolution bits (step 72). An example is         given below, from U.S. patent application Ser. No. 11/754,965,         which is entitled “Quantization of Channel State Information in         Multiple Antenna Systems”.     -   4. Store the description of the Voronoi regions and centroids         for m-tier of the VQ (step 74).     -   5. If m is smaller than M (step 76), continue the design in the         following way:         -   a) Create a large list of possible channel realizations             (step 78).         -   b) Using the m-tier VQ, quantize (step 80) the above channel             realizations.         -   c) Select channel realizations corresponding to each of the             m-tier VQ indices (step 82).         -   d) Within each of the m-tier Voronoi regions, perform m+1             tier design of a vector quantizer using any of various VQ             designs, such as in U.S. patent application Ser. No.             11/754,965 and shown below. The algorithm uses N_(m+1)             resolution bits within each region and forces one of the             tier m+1 centroids within each region to be identical to the             m-tier centroid corresponding to this region (step 84). The             codebook entries of the new codebook for each region are now             considered to be the group of tier m+1 entries associated             with the tier m codebook entry for that region.         -   e) The reused m-tier centroid is assigned N_(m+1) index bits             equal to 0.         -   f) Increase m by 1 (step 86).         -   g) Go to step 4.     -   6. If m≧M, finish the VQ design (step 76).     -   7. Design the modulation matrices corresponding to the designed         multi-tiered VQ using any of various modulation matrix design         techniques such as the algorithm shown in U.S. patent         application Ser. No. 11/754,965, as shown below (step 88). The         algorithm is now done (step 90).

The rationale behind re-using one of the m tier centroids at design phase of (m+1)-tier centroids and Voronoi regions (see bullet 5d above) is that the same set of modulation matrices can be used in a system where different users report their quantized channel information using different VQ tiers. As all m-tier centroids are contained in (m+1)-tier centroids, the effective indices can be easily used to decide which centroid must be used. FIG. 6 shows an example: a single primary region with index 1111 92 and secondary regions with indices of 000, 001, etc 94 thus giving them effective indices of 1111000, 1111001, etc. 96 Note that the primary centroid is in the same place as the secondary centroid with index 000. Moreover, thanks to such a design, the base station may support users with different implementations of the vector quantizers, e.g., varying number of VQ tiers. Thanks to the embedding of the codewords, all such situations will be supported.

The algorithm from “Quantization of channel state information in multiple antenna systems” is as follows:

For the case of a single receiver active at a time, we introduce a heuristic distortion metric which is expressed as

γ_(V)(n;H)=∥DV ^(H) {circumflex over (V)}(n)−D∥ _(F)  (1)

where {circumflex over (V)}(n) is the nth entry in the predefined set of channel diagonalization matrices and ∥·∥_(F) is the Frobenius norm. We omitted subscript entries j in (1) for the clarity of presentation.

We assume that n=0, 1, . . . 2^(N) ^(V) −1 where N_(V) is the number of bits per channel realization in the feedback link needed to represent the vectors {circumflex over (V)}(n). To design the quantizer using (1), we divide the whole space of channel realizations H into 2^(N) ^(V) regions V_(i) where

V _(i) ={H:γ _(V)(i;H)<γ_(V)(j;H) for all j≠i}.  (2)

The algorithm starts by creating a codebook of centroids {circumflex over (V)} and, based on these results divides the quantization space into regions V_(i). The codebook is created as follows:

-   -   1. Create a large training set of L random matrices H(l).     -   2. For each random matrix H(l), perform singular value         decomposition to obtain D(l) and V(l) as

y _(j) =H _(j) x _(j) +n _(j)=(U _(j) D _(j) V _(j) ^(H))(V _(j) x _(j))+n _(j)  (3)

-   -   3. Set iteration counter i=0. Create a set of 2^(N) ^(V) random         matrices Ĥ(n).     -   4. For each matrix Ĥ(n) calculate corresponding {circumflex over         (V)}^((i))(n) using singular value decomposition.     -   5. For each training element H(l) and codebook entry {circumflex         over (V)}^((i))(n) calculate the metric in (1). For every l         choose indexes n_(opt)(l) corresponding to the lowest values of         γ_(V)(n;H(l)).     -   6. Calculate a new set {circumflex over (V)}^((i+1))(n) as a         form of spherical average of all entries V(l) corresponding to         the same index n using the following method. (The direct         averaging is impossible since it does not preserve orthogonality         between eigenvectors.) For all n calculate the subsets         L(n)={l:n_(opt)(l)=n} and if their respective cardinalities         |L(n)|≠0 the corresponding matrices Q ^((i+1))(n) can be         obtained as

$\begin{matrix} {{{\overset{\_}{Q}}^{({i + 1})}(n)} = {\frac{1}{{L(n)}}{\sum\limits_{l \in {L{(n)}}}{{v(l)}^{1}{{OV}(l)}^{H}}}}} & (4) \end{matrix}$

-   -    where ¹O is an n_(T)×n_(T) all-zero matrix with the exception         of the upper-left corner element equal to 1. Finally, using         singular value decomposition, calculate {circumflex over         (V)}^((i+1))(n) from

Q ^((i+1))(n)={circumflex over (V)} ^((i+1))(n)W({circumflex over (V)} ^((i+1))(n))^(H)  (5)

-   -    where W is a dummy variable.     -   7. Calculate the average distortion metric γ _(V)         ^((i+1))=1/LΣ_(l)γ_(V)(n_(opt)(l);H(l)).     -   8. If distortion metric fulfils | γ _(V) ^((i+1))− γ _(V)         ^((i))|/ γ _(V) ^((i))<Θ, stop. Otherwise increase i by 1 and go         to 5).

Upon completion of the above, algorithm, the set of vectors {circumflex over (V)} can be used to calculate the regions in (2).

Having optimized power-independent entries in the codebook of channel eigenmode matrices {circumflex over (V)}, the next step is to create a codebook for power allocation Ŝ. We use a distortion metric defined as

$\begin{matrix} {{{\gamma \;}_{S}\left( {k;H;P} \right)} = \begin{matrix} {\det \left\lbrack { + {HQH}^{H}} \right\rbrack} \\ {\det\left\lbrack { + {H\; {\hat{V}\left( n_{opt} \right)}{\hat{S}(k)}{{\hat{V}}^{H}\left( n_{opt} \right)}H^{H}}} \right\rbrack} \end{matrix}} & (6) \end{matrix}$

where Ŝ(k) is the kth entry in the predefined set of channel water-filling matrices and {circumflex over (V)}(n_(opt)) is the entry in {circumflex over (V)} codebook that minimizes metric (1) for the given H. We use k=0, 1, . . . 2^(N) ^(S) −1 where N_(S) is the number of bits per channel realization in the feedback link needed to represent the vectors Ŝ(k). Minimizing the metric in (6) is equivalent to minimizing the capacity loss between the optimum water-filling using Q and the quantized water-filling using {circumflex over (V)} and Ŝ.

Similarly to the previous problem, we divide the whole space of channel realizations H into 2^(N) ^(S) regions S_(i)(P) where

S _(i)(P)={H:γ _(S)(i;H;P)<γ_(S)(j;H;P) for all j≠i}.  (7)

and to create the codebook Ŝ, we use the following method:

-   -   1. Create a large training set of L random matrices H(l).     -   2. For each random matrix H(l), perform water-filling operation         to obtain optimum covariance matrices Q(l) and S(l).     -   3. Set iteration counter i=0. Create a set of 2^(N) ^(S) random         diagonal matrices Ŝ^((i))(k) with Tr(Ŝ^((i))(k)=P.     -   4. For every codebook entry Ŝ^((i))(k) and matrix Q(l),         calculate the metric as in (6). Choose indexes k_(opt)(l)         corresponding to the lowest values of γ_(s)(k;H(l);P).     -   5. If γ_(S)(k_(opt)(l); H(l);P)>γ_(eq)(H(l);P) where         γ_(eq)(H(l);P) is the metric corresponding to equal-power         distribution defined as

$\begin{matrix} {{\gamma_{eq}\left( {{H(l)};P} \right)} = \begin{matrix} {\det \left\lbrack { + {{H(l)}{Q(l)}{H^{H}(l)}}} \right\rbrack} \\ {{\det \left\lbrack { + {{P/n_{T}}{H(l)}{H^{H}(l)}}} \right\rbrack},} \end{matrix}} & (8) \end{matrix}$

-   -    set the corresponding entry k_(opt)(l)=2^(N) ^(S) . For all k         calculate the subsets L(k)={l:k_(opt)(l)=k}.     -   6. For all k=0, 1, . . . 2^(N) ^(S) −1 for which |L(k)|≠0,         calculate a new set Ŝ^((i+1))(k) as the arithmetic average

$\begin{matrix} {{{\hat{S}}^{({i + 1})}(k)} = {\begin{matrix} 1 \\ {{L(k)}} \end{matrix}{\sum\limits_{l \in {L{(k)}}}{{S(l)}.}}}} & (9) \end{matrix}$

-   -   7. Calculate the average distortion metric

$\begin{matrix} {{\overset{\_}{\gamma}}_{S}^{({i + 1})} = {\frac{1}{L}{\sum\limits_{l}{\min {\left\{ {{\gamma_{S}\left( {{k_{opt}(l)};{H(l)};P} \right)},{\gamma_{eq}\left( {{H(l)};P} \right)}} \right\}.}}}}} & (10) \end{matrix}$

-   -   8. If distortion metric | γ _(S) ^((i+1))− γ _(S) ^((i))|/ γ         _(S) ^((i))<Θ, stop. Otherwise increase i with 1 and go to 4).

The set of vectors Ŝ is then used to calculate the regions in (7). Since water-filling strongly depends, on the power level P and {circumflex over (V)}, optimally the Ŝ should be created for every power level and number of bits N_(V) in eigenvector matrix codebook. (2).

In the multi-user case, we follow the approach of Spencer et al, where each user performs singular value decomposition of H_(k)=U_(k)S_(k)V_(k) ^(H) and converts its respective H_(k) to a n_(T)-dimensional vector h_(k) as

h _(k) =u _(k) ^(H) H _(k) =s _(k) ^(max) v _(k) ^(H)  (11)

where s_(k) ^(max) is the largest singular value of S_(k) and u_(k) and v_(k) are its corresponding vectors from the unitary matrices U_(k) and V_(k), respectively.

We use the linear block diagonalization approach which eliminates. MUI by composing the modulation matrix B[S] of properly chosen null-space eigenmodes for each set S. For each receiver iεS, the ith row of the matrix. H[S] is first deleted to form H[S_(i)]. In the next step, the singular value decomposition is performed to yield H[S_(i)]=U[S_(i)]S[S_(i)]V^(H)[S_(i)]. By setting the ith column of B[S] to be equal to the rightmost vector of V[S_(i)] we force the signal to the ith receiver to be transmitted in the null-space of the other users and no MUI will appear. In other words, the channel will be diagonalized with d_(i), being the entries on the diagonal of H[8]B[8]. This leads to formula

$\begin{matrix} {R^{linear} = {\max\limits_{S}{\sum\limits_{i = 1}^{n_{\Gamma}}\left\lbrack {\log_{2}\left( {{\xi \lbrack S\rbrack}d_{i}^{2}} \right)} \right\rbrack_{+}}}} & (12) \end{matrix}$

where ξ[S] is the solution of the water-filling equation.

We assume that N_(μ) is the number of bits per than channel realization in the feedback link needed to represent the vectors v_(k) in (11). We divide the space of all possible v's into 2^(N) ^(μ) regions μ_(i)

μ_(i) ={v:γ _(μ)(i;v)<γ_(μ)(j;v) for all j≠i}  (13)

where γ_(μ)(n;v) is a distortion function. Within each region μ_(i) we define a centroid vector {circumflex over (v)}(i), which will be used as a representation of the region. The design of the codebook {circumflex over (v)} can be done analytically and/or heuristically using for example the Lloyd algorithm. In this work, we define the distortion function as the angle between the actual vector v and {circumflex over (v)}(i): γ_(μ)(i;v)=cos⁻¹({circumflex over (v)}(i)·v), which has been shown by Roh and Rhao to maximize ergodic capacity, and use Lloyd algorithm to train the vector quantizer. Note that the construction of {circumflex over (v)} is independent of the transmit power.

We assume that N_(s) is the number of bits per channel realization in the feedback link needed to represent the scalar s_(k) ^(max) in (11). We divide the space of all possible channel realizations s=s^(max) into 2^(N) ^(s) regions s_(i)

s _(i) ={s:|ŝ(i)−s|<|ŝ(j)−s| for all j≠i}  (14)

where ŝ(i) are scalar centroids representing regions s_(i). In this work, we perform the design of the codebook ŝ using the classical non-uniform quantizer design algorithm with distortion function given by quadratic function of the quantization error as ε(i;s)=(s−ŝ(i))².

The construction of the codebook ŝ is generally dependent on the transmit power level. However, the differences between the codebooks ŝ for different power regions are quite small. This allows us to create only one codebook ŝ and use it tar all transmit powers.

The calculation of the modulation matrix {circumflex over (B)} is based on the given codebook {circumflex over (v)}. We assume that the quantization of the channel eigenmodes is performed at the receiver side and each user transmits back its codebook index i_(k). The indices are then used at the transmitter side to select the modulation matrix {circumflex over (B)}(i₁, i₂, . . . i_(K)). Since, from the linear transmitter point of view, ordering of the users is not important, we use the convention that the indices (i₁, i₂, . . . i_(K)) are always presented in the ascending order. For example, in a system with K=2, n_(T)=2 and 1-bit vector quantizers {circumflex over (v)}, there will exist only three possible modulation matrices corresponding to sets of {circumflex over (v)} indices (1, 1), (1, 2) and (2, 2).

In the context of vector quantizing, the design of the modulation matrices can no longer be based on the algorithm presented for the single user case. Using this method with quantized versions of h_(k) produces wrong result when identical indices i_(k) are returned and the receiver attempts to jointly optimize transmission to the users with seemingly identical channel vectors ĥ_(k). Instead, we propose the following algorithm to optimize the set of matrices {circumflex over (B)}(i₁, i₂, . . . i_(K)):

-   -   1. Create a large set of Nn_(T) random matrices H_(k), where N         is the number of training, sets with n_(T) users each.     -   2. For each random matrix H_(k), perform singular value         decomposition and obtain h_(k) as in (11).     -   3. For each vector h_(k) store the index i_(k) of the         corresponding entry {circumflex over (v)}(i_(k)).     -   4. Divide the entire set of matrices H_(k) into N sets with         n_(T) elements each.     -   5. Sort the indices i_(k) within each set l in the ascending         order. Map all unique sets of sorted indices to a set of unique         indices I_(B) (for example (1,1)→I_(B)=1; (1,2)→I_(B)=2;         (2,2)→I_(B)=3 . . . ).     -   6. In each set l, reorder the corresponding channel vectors         h_(k) according to their indices i_(k) and calculate the optimum         B_(l) using the block diagonalization method described above.     -   7. Calculate a set {circumflex over (B)}(I_(B)) as a column-wise         spherical average of all entries B_(l) corresponding to the same         index I_(B).

After calculation of |I_(B)| modulation matrices {circumflex over (B)}, the remaining part of system design is the calculation of the water-filling matrices {circumflex over (D)}, which divide the powers between the eigenmodes at the transmitter. The procedure for creation of codebook {circumflex over (D)} is similar to the above algorithm, with the difference that the entries ŝ(n_(k)) are used instead of {circumflex over (v)}(i_(k)), and the spherical averaging of the water-filling matrices is performed diagonally, not column-wise. Explicitly:

-   -   1. Create a large set of Nn_(T) random matrices H_(k), where N         is the number of training sets with n_(T) users each.     -   2. For each random matrix H_(k), perform singular value         decomposition and obtain h_(k) as in (11).     -   3. For each vector h_(k) store the index n_(k) of the         corresponding entry ŝ(n_(k)).     -   4. Divide the entire set of matrices H_(k) into N sets with         n_(T) elements each.     -   5. Sort the indices n_(k) within each set l in the ascending         order. Map all unique sets of sorted indices to a set of unique         indices I_(D) (for example (1,1)→I_(D)=1; (1,2)→I_(D)=2;         (2,2)→I_(D)=3 . . . ).     -   6. In each set l, reorder the corresponding channel vectors         h_(k) according to their indices n_(k) and calculate the optimum         D_(l) using the method of waterfilling of (12).     -   7. Calculate a set {circumflex over (D)}(I_(D)) as a diagonal         spherical average of all entries D_(l) corresponding to the same         index I_(D).

Referring to FIGS. 7 and 8, based on the design of the multi-tiered codebooks D and V as in the previous section, the system will operate as follows:

-   -   1. Initialize transmission epoch to t=1.     -   2. Set m_(k)=1 at each receiver k, (all users will use separate         indices m_(k)) (step 112).     -   3. Set m^(k)=1 separately for each receiver at the transmitter         side. The transmitter-side indices m^(k) should be mapped to         their respective receiver-side indices m_(k) (step 110).     -   4. Each receiver estimates its channel matrix H[t] (step 40).     -   5. Each receiver performs the vector quantization of the channel         using m=1 tier quantizers described above (step 114).     -   6. The m-tier N_(m)-bit long indices are fed back to the         transmitter.     -   7. The transmitter performs the selection of active users using         any method (maximum fairness, maximum throughput etc.) and         chooses the optimum modulation matrices using a VQ method such         as the method described in U.S. patent application Ser. No.         11/754,965.     -   8. The signal is transmitted to the selected active receivers.     -   9. Increase transmission epoch as t=t+1.     -   10. Each receiver estimates its channel matrix H[t] (step 40).     -   11. Each receiver performs the vector quantization of the         channel using m-tier quantizers described above (step 114).     -   12. Each receiver that recognizes (step 116) that its quantized         channel's m_(k)-tier Voronoi region in the t+1 epoch is         identical to the m_(k)-tier Voronoi region in epoch t performs         the following steps:         -   a) Unless 118 m_(k)=M, increase the receiver's index to             m_(k)=m_(k)+1 (step 120).         -   b) The channel realization within the unchanged Voronoi             region is quantized using the new m_(k) tier quantizer (step             114).         -   c) The receiver uses a known mechanism (see later in the             document) to signal to the transmitter the new m_(k)-tier of             VQ.         -   d) The N_(mk) bits long indices are fed back to the             transmitter (step 124).         -   e) Transmitter increases its index as m^(k)=m^(k)+1 (step             126).     -   13. Each receiver that recognizes (step 116) that its quantized         channel's m_(k)-tier Voronoi region in the t+1 epoch is not         identical to the m_(k)-tier Voronoi region in epoch t performs         the following steps:         -   a) Unless (step 128) m_(k)=1, decrease the receiver's index             m_(k) to the last tier where the Voronoi regions m_(k-1) are             identical in both t and t+1 epochs (step 130).         -   b) If no such tier can be found, set m_(k)=1, otherwise             update the m_(k) to a value for which Voronoi regions             m_(k−1) are identical (step 130).         -   c) The channel realization is quantized using the m_(k)-tier             quantizer (step 114).         -   d) The receiver uses a known mechanism (see below) to signal             to the transmitter the new m_(k) tier of VQ.         -   e) The N_(mk) bits long indices are fed back to the             transmitter (step 124).         -   f) Transmitter decreases its index as m^(k)=m_(k) (step             134).     -   14. The transmitter selects the modulation matrices based on the         indices fed from the receivers and each receiver's separate         m^(k) index stored at the transmitter side (step 136).     -   15. The modulation matrices are used to transmit the information         to the selected receivers (step 138).         The example of the algorithm's operation shown in FIGS. 7-8 for         one mobile receiver uses M=3 tiered quantizer. In this scenario,         CSI vector stays in tier-1 Voronoi region 20 in first 5 frames         F1-F5, in the tier-2 region 24 in first 5 frames F1-F5, and         tier-3 region 100 in frames F2 and F3. Receiver R1 recognizes         the subsequent Voronoi regions and adjusts the tier (m) 104 of         the used quantizer accordingly by increasing and decreasing         quantizer resolution. The quantizer indices each representing         the centroid 22,26,102 of its respective Voronoi region in the         appropriate tier, are then fed to the base station B1 that         combines them properly so that the effective CSI resolution N         varies in time (t increases from F1-F6) depending on the rate of         channel changes. Base station B1 chooses modulation matrix 56         and transmits signal 60. At each time the resolution is equal to         that of an untiered quantizer with a number of bits equal to the         number of bits 106 representing the tier used (N_(m)) plus that         for all lower tiers. It can be clearly seen that the proposed         algorithm allows the system to automatically adjust the         resolution to the speed of channel changes.         FIG. 9 shows a graphical representation of the algorithm.

The set of indices m_(k) at the receivers should be matched to the indices m^(k) at the transmitter. If the transmitter uses the index m^(k) that corresponds to the wrong m_(k)-tier of the receiver VQ, the resulting loss of performance may be very significant. In general, the index of the quantized channel vector at the transmitter is reconstructed as:

AAAABBBCCC . . . .

where AAAA corresponds to m=1 tier N₁ indexing bits, BBB corresponds to m=2 tier N₂ indexing bits etc. (see FIG. 6 for an example). At any given time, the transmitter receives only m-tier index bits (AAAA, BBB, CCC etc.) and it must be able to establish which tier those bits correspond to. For example, there must be a signaling method allowing the transmitter to distinguish between two consecutive transmissions such as BBB, BBB where the channel vector moved away from one tier-2 centroid to another, from the BBB, CCC transmission, where the channel vector stayed in the same tier-2 region BBB and tier-3 quantization was used in the CCC word.

Various methods may be used for transmitting index information from the receiver to the transmitter such as:

-   -   1. Direct indexing of VQ words. In order to let the transmitter         know, which VQ tier is used, the actual VQ codeword index is         extended with the bit representation of the index m_(k) of each         receiver. Example:     -   MMAAAA, MMBBB, MMCCC . . . .     -   Where, for example, two bits MM are used to represent one of the         four quantization tiers in the system. The drawback of this         method is that the additional feedback load is required to         transmit information about m_(k) indices.     -   2. Varying length of different m tier VQ words. In this method,         each m-tier of the CSI quantizer is characterized by different         number of indexing bits N_(m). Example:     -   AAAA, BBB, CC, D . . . .     -   where four bits AAA are used to represent tier-1 quantization, 3         bits BBB are used to represent tier-2 quantization etc. The         advantage of this system is that there is no need to transmit         additional bits M as in the previous method. The drawback of         this method is that there must be another mechanism allowing the         transmitter to count how many bits were actually sent from the         receiver and the varying feedback load.     -   3. Channel prediction based assessment of VQ tier. In this         method, each m-tier of the CSI quantizer may be characterized by         any number of bits N_(m) and the statistical channel         characterization is used by the transmitter to decide whether         the channel vector stayed in the previous m-tier Voronoi region         or moved away from it. The advantage of this system is that it         leaves a large degree of freedom for designing the feedback         link. The drawback of this method is that the complexity of         transmitter design grows and there may be erroneous decisions on         the tier of quantizer used by the receivers.     -   4. Hybrid solutions combining the previous three methods in any         way that is suitable from system design point of view.

In the course of the system operation, it may happen that some of the transmitter indices M^(k) will no longer be synchronized with corresponding receiver indices m_(k). Such a situation will typically happen when one of the feedback messages from a receiver has not been detected at the transmitter (i.e., the transmitter lacks channel quantization index for the current transmission epoch) or the received message with the indexing information does not agree with the expected quantization tier m.

In practical communication systems, two erroneous situations can occur:

-   -   The received feedback message shows the situation when         m_(k)>m^(k)+1. Such a situation is not allowed during the course         of the normal operation since the receiver may only step back to         the lower tier quantizers or increase the current one by 1.     -   The transmitter did not receive any feedback information due to         the feedback link problems.

Various methods may be used to solve the problem such as:

-   -   1. Use channel prediction to recover the incorrect index. The         previously used indices are used to extrapolate the actual         channel information index.     -   2. Deactivate the user for the next transmission epoch and send         the VQ RESET message to it. If the transmitter cannot reliably         decide, which channel index was reported by the receiver, it         sends a special VQ RESET message to the receiver containing the         last value of the effective index at the transmitter AAAABBBCC .         . . without the last tier of bits (in other words, bits up to         the level m^(k)−1 are communicated to the base station). The         receiver than establishes, whether the same tier of the VQ can         be used or whether it has to step back to a lower tier. The new         indices are sent to the transmitter and the system resumes the         usual operation.

FIG. 10 shows a typical set of curves representing the eigenmode and singular value coherence times for a 2×2 MIMO system. For definitions see B. Mielczarek and W. Krzymien, “Influence of CSI feedback delay on capacity of linear multi-user MIMO systems,” in Proc. IEEE WCNC., Hong Kong, March 2007, pp. 1188-1192. As one can see, the length of time in which the first tier Voronoi regions do not change decreases with increasing resolution and normalized Doppler frequency of the channel f_(D)T_(frame). For example, with f_(D)T_(frame)=0.02, the Voronoi regions of the eigenmode quantizers will stay the same for approximately 6 consecutive frames when the eigenmode quantizer uses N=4 bit resolution. When the quantizer uses N=7 bits, only 2-3 consecutive frames will have identical first tier Voronoi regions. In general, in order to improve system's throughput, it is preferable to use higher resolution of VQ but the price for the improvement is the high feedback burden and frequent changes of the high resolution indices. As shown in example below, by using a multi-tier VQ design, it is actually possible to achieve very good system performance and significantly reduce the required feedback bit rate.

In FIG. 11, the simulation results are shown for a system with 2-tier eigenmode quantizer (M=2) using N₁=4 and N₂=3. The system has been designed using algorithms described in this application and starts by transmitting 4 bits (AAAA) and, if the 1-tier Voronoi regions are the same in the consecutive frames, only 3 bits (BBB) are sent to the transmitter. We compare them with 1-tier systems with N₁=7 and N₁=4. The implemented algorithms are tested on a system with 2 base station antennas and 10 users with 2 receive antennas each. We test three channels with maximum normalized Doppler frequencies equal to 0.01, 0.02 and 0.1 as in FIG. 11. As one can see, the throughput gap between the conventional 1-tier vector quantizers for N₁=7 and N₁=4 is quite large (around 3 dB at 10 bpcu)—any increase of throughput in such simple systems requires increase of feedback bandwidth. However, if the channel is assumed to have memory, by using our proposed approach, it is possible to attain almost the same performance with multi-tier CSI quantization. In our example, the maximum feedback burden is set to 4 bits/frame/receiver for the 2-tier system but the performance is almost equivalent to 7 bit feedback system for a large range of Doppler frequencies. Hence, by proper choice of the number of tiers and their corresponding resolutions, it is possible to design the practical systems for a wide variety of channel conditions, required throughput performance and maximum feedback link bit rates.

In the claims, the word “comprising” is used in its inclusive sense and does not exclude other elements being present. The indefinite article “a” before a claim feature does not exclude more than one of the feature being present. Each one of the individual features described here may be used in one or more embodiments and is not, by virtue only of being described here, to be construed as essential to all embodiments as defined by the claims.

Immaterial modifications may be made to the embodiments described here without departing from what is covered by the claims. 

1. A method of quantizing channel state information in a multiple-input transmission system having at least a transmitter and a receiver, the method comprising the steps of: quantizing information concerning a first channel state to produce first quantized information, and sending the first quantized information to the transmitter; and quantizing information concerning a second channel state by reference to (1) the second channel state and (2) the first quantized information, to produce a second quantized information, and sending the second quantized information to the transmitter. 2-20. (canceled)
 21. An apparatus, comprising: a receiver comprising one or more antennas and a processing module processing module, said module configured to carry out the method steps of claim 1; and a transmitter comprising antennas and a processing module configured to process channel state information quantized by the method of claim 1 from one or more receivers and modulate a signal sent to one or more of the receivers according to the channel state information. 